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£ . ABSTRACT 

cd 

, We analyse the iGband Hubble diagram for a sample of brightest cluster galax- 

00 ■ ies (BCGs) in the redshift range < z < 1. In good agreement with earlier stud- 

Cn| \ ies, we confirm that the scatter in the absolute magnitudes of the galaxies is small 

(0.3 magnitudes). The BCGs exhibit very little luminosity evolution in this redshift 
range: if qo — 0.0 we detect no luminosity evolution; for qo = 0.5 we measure a small 
negative evolution (i.e., BCGs were about 0.5 magnitudes fainter at z = 1 than to- 
day). If the mass in stars of these galaxies had remained constant over this period 
of time, substantial positive luminosity evolution would be expected: BCGs should 
have been brighter in the past since their stars were younger. A likely explanation 
f- ^ for the observed zero or negative evolution is that the stellar mass of the BCGs has 

OO I been assembled over time through merging and accretion, as expected in hierarchi- 

Q\ . cal models of galaxy formation. The colour evolution of the BCGs is consistent with 

that of an old stellar population (zf orrn > 2) that is evolving passively. We can thus 
use evolutionary population synthesis models to estimate the rate of growth in stellar 
mass for these systems. We find that the stellar mass in a typical BCG has grown by 
a factor ~ 2 since z ~ 1 if go = 0-0 or by factor ~ 4 if go = 0.5. These results are 
, in good agreement with the predictions of semi-analytic models of galaxy formation 

and evolution set in the context of a hierarchical scenario for structure formation. The 
models predict a scatter in the luminosities of the BCGs that is somewhat larger than 
, the observed one, but that depends on the criterion used to select the model clusters. 

•i-H . 
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1 INTRODUCTION 



Brightest cluster galaxies (BC Gs) have be en exte nsively 



studied at optical wavelengths ((Peach 1970), 1972; (Glunn 


& Okc 1975), (Sandagc, Kristian & Westphal 1976); ([Kris- 


tian, Sandage 


& Westphal 19780; (IKristian, Sandage & 


Westphal 1978 


); (Hoessel 1980f); ([Schneider, Gunn & Hoes- 


scl 1983a ),b; (| 


jauer & Postman 1992|); (postman & Lauer 


1995)). The small scatter in their absolute magnitudes (< 



0.3 magnitudes up to z ~ 0.5, see e.g. ( pandage 19 
their high luminosities make them useful standard candles 
for classical cosmological tests involving the Hubble dia- 
gram, such as the determination of qo and the variation of 
Ho with redshift. However, there is now firm evidence that 
significant evolution has taken place in the colours and opti- 
cal luminosities of early-ty pe galaxies (including the BCGs ) 
from z ~ to z ~ 1 (e.g. ( Aragon-Salamanca et al. 1993 ); 



(Rakos fc Schombert 1995); (Oke, Gunn fc Hoessel 1996) 



(Lubin 1996); (Stanford, Eisenhardt & Dickinson 1997)) im- 
plying that the Hubble diagram could be seriously affected 
by evolutionary changes even at relatively modest redshifts. 
Conclusions concerning go cannot be derived from such dia- 
grams until these changes are well-understood. Indeed, the 
Hubble diagram for the BCGs may well have more to teach 
us about galaxy evolution than about cosmology. 

The study of the BCGs in the near-infrared K-band 
(2.2fi) has two main advantages over optical studies: first, 
the fc-corrections are appreciably smaller (indeed, they are 
negative) and virtually in dependent of the spectral typ e 
of the galaxies (see, e.g., (Aragon-Salamanca et al. 1993), 
1994); and second, the light at long wavelengths is domi- 
nated by long-lived stars. Thus the JT-band luminosity is a 
good measure of the total stellar mass in the galaxies. 

The A"-band Hubble diagram has long been used as an 
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evolutionary dia gn ostic for radio galaxies ((|Grasdalen 1980|); 
flLebofsky 198C| ); ( |Lebofsky fc Eisenhardt 1986| ); dLiiiy& 
Longaiil 1982|), 1984~7|Lilly 1989a|),b), which also have a rel- 



atively mall intrinsic dispersion in their luminosities (~ 0.4 



magnitudes for 1< z < 2.2; ( [Lilly 19894 ). Substantial pos- 
itive luminosity evolution was found: radio galaxies were 
about one magnitude brighter at z ~ 1 than they are to- 
day. Aragon-Salamanca et al. (1993) measured _K"-band lu- 
minosities for a sample of 19 BCGs in the redshift range 
< z < 1 and found that the K magnitude-redshift relation 
is tighter (scatter = 0.3 magnitudes) than for radio galaxies. 
Moreover, these authors did not detect significant K lumi- 
nosity evolution for the BCGs, and concluded that their evo- 
lutionary properties were different from those of radio galax- 
ies: powerful radio galaxies are apparently less homogeneous 
in their evolutionary behaviour and show stronger luminos- 
ity evolution. Aragon-Salamanca et al. also found that the 
colours of the early-type galaxies (including the BCGs) in 
these clusters have evolved since z ~ 1 at a rate consistent 
with a scenario in which the bulk of their stars formed at 
relatively high redshifts (z > 2) and evolved passively there- 
after. This conc lusion has been confirmed b y several more 
recent studies (( Etakos fc Schombert 1995); (Dke, Gunn & 
HoesselT T996|); ([Lubin 1996|); QEUis et al. 1997| ); ^Stanford 



rdt fc Dickinson 1997)) 



Eisenh; 

The recent development of semi-analytic techniques has 
provided theorists with the tools to make predictions for the 
formation and evolution of galaxies, using physically moti- 
vated models set in the contex t of hie rarchical structure for- 
mation in the universe (e. g., (|White 

199l |)~|Fac"ey fc Silk 199l|); ([Kauffmann, White fc Guider- 



( |Kauffmann7Guiderdoni fc White 1994| ); QCole 
^eyT*et**^iL*T99£ )). The properties of galaxies 



doni 193, 
et al. 1 )9 

in these models are in broad agreement with the present day 
characteristics of galaxies, such as the distribution of lumi- 
nosities, colours and morphologies and with the properties of 
galaxies at high redshift, including the f aint galaxy counts , 
colours and redshift distributions (e.g. (Kauffmann 1995); 



(Kauffmann 1996); (Kauffmann fc Chariot 1997); (Baugh 



Cole fc |Frenk 1996a| ); (|Baugh et. al. 1997|) ). In this paper, 
we will test whether the models also predict the right be- 
haviour for a special class of galaxies, namely the BCGs. In 
the models, these galaxies grow in mass both from the cool- 
ing of gas from the surrounding hot halo medium and from 
the accretion of "satellite" galaxies that fall to the cluster 
centre as result of dynamical friction and then merge. 

In section 2 we present the observational results and de- 
rive the luminosity evolution of the BCGs in the sample. In 
section 3 we briefly describe two independent semi-analytic 
models of galaxy formation, outline their basic assumptions, 
limitations and uncertainties, and discuss the predictions 
they make for the evolution of the BCGs. Section 4 com- 
pares the model predictions with the observational results. 



2 OBSERVATIONAL RESULTS 
2.1 Photometric data 

The clusters in our sample are all optically-selected, al- 
though most of them have also been detected in X-rays. At 
z < 0.37 they come from the Abell (1957) and Abell, Corwin 



& Olowin (1989) catalogues. At higher redshifts they have 
been found from the projected density contrast in deep opti- 
cal imag ing surveys ((Gunn, Hoessel & Oke 1986); (Couch et 



al. 1991)) and should represent the richest clusters present 



at each redshift. 

The A'-band data analysed here come from the photo- 
metric study carried out by Aragon-Salamanca et al. (1993) 
of galaxy clusters in the < z < 1 redshift range (19 
BCGs) complem ented with new _K"- band images of 3 clus- 
ters at z ~ 3 ((Barger et al. 1996)) and 6 clusters with 



0.39 < z < 0.56 ( QBarger et al. 1997[ )). The images were ob- 
tained at the 3.8-m UK Infrared Telescope, the 3.9-m Anglo- 
Australian Telescope and the Palomar Observatory 1.52-m 
Telescope. There is some overlap in the clusters observed 
by Aragon-Salamanca et al. and Barger et al., bringing the 
number of BCGs studied here to a total of 25 (see Table 0). 
For the clusters in common, both sets of photometry agree 
very well (within the estimated errors) and we have chosen 
the better quality image in each case. 

We refer the reader to the original papers for a detailed 
description of the data. We will concentrate here on the 
aspects relevant to this paper. The photometric zero-point 
uncertainties in the if-band images were always small (0.02- 
0.04 magnitudes) and are completely negligible for our anal- 
ysis. As in Aragon-Salamanca et al. (1993), K photometry 
for the BCGs was obtained inside a fixed metric aperture 
of 50 kpc diameter^] so that the minimum aperture at high 
redshift is ~ 5 arcsec thus minimising the effect of seeing. 
Because of the ambiguity of what constitutes a "galaxy" 
(~ 30% of local first-ranked cluster galaxies are multiple- 
nucleus systems: i.e., more than one "ga laxy" occurs within 
our metric aperture — see (Hoessel 198C)) we follow the ap- 
proach of Schneider, Gunn & Hoessel (1983a) and adopt a 
working definition of the brightest cluster galaxy as the re- 
gion of maximum cluster light enclosed in our metric aper- 
ture. The photometry data from the AAT and PO was trans- 
formed into the K-banA system of the UKIRT (where most 
of the data comes from) using colour terms determined from 
the known filter responses. The colour terms were always 
< 0.08 magnitudes, and we estimate that the uncertainty 
in the transformation is well below 2% and thus negligible. 
Galactic reddening corrections were estimated from Burstein 
& Heiles' (1982) maps, but since the galactic latitudes were 
always reasonably high, they were very small in the if-band. 
Finally, small seeing corrections and fc-corrections were also 
applied as in Aragon-Salamanca et al. Table |l| contains the 
corrected if-magnitudes for the BCG sample. Values for 
go = 0.0 and 0.5 are given to account for the difference in 
projected aperture at a given metric aperture as a function 
of the deceleration parameter. We estimate that the total 
uncertainty of the corrected magnitudes is below 10%, and 
will contribute very little to the observed scatter (see original 
data papers for a detailed discussion of the uncertainties). 

Figure |l| shows the K magnitude-redshift relation 
(Hubble diagram) for the BCGs in our sample. The solid 
lines show no-evolution predictions which only take into ac- 
count the distance modulus as a function of z, normalised 



* Ho = 50 kms" 1 Mpc~ 1 assumed throughout. Because of the 
uncertainty in the value of go we w iU carry out parallel analyses 
for go = 0.0 and 0.5 
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Table 1. Corrected rest-frame _R"-band magnitudes for the BCGs 



Cluster 


z 


K° 

"q =0-0 




reference 


Abell 1656 a (NGC4889) 


0.023 


8.89 


8.88 


(1) 


Abell 2199 (NGC6166) 


0.030 


9.50 


9.49 


(1) 


Abell 2197 (NGC6173) 


0.031 


9.54 


9.53 


(1) 


Abell 2151 (NGC6034) 


0.037 


10.49 


10.38 


(1) 


A 1 11 nno 

Abell 963 


0.206 


13.61 


13.55 


(1) 


Abell 1942 


0.224 


13.92 


13.88 


(1) 


AC103 


0.310 


14.62 


14.56 


(2) 


AC114 


0.310 


14.54 


14.46 


(2) 


AC118 


0.310 


14.77 


14.67 


(2) 


C12244— 02 


0.329 


15.35 


15.31 


(1) 


A 1 ii o 

Abell 370 


0.374 


15.05 


14.95 


(1) 


C10024+16 


0.391 


15.29 


15.14 


(1).(3) 


3C295 


0.460 


15.26 


15.14 


(3) 


C10412— 65 b 


0.510 


15.60 


15.55 


(3) 


C11601+42 


0.540 


16.23 


16.13 


(3) 


C10016+16 


0.546 


16.21 


16.02 


(1),(3) 


C10054-27 c 


0.563 


16.50 


16.37 


(1),(3) 


C10317+1521 


0.583 


17.08 


17.01 


(1) 


C10844+18 d 


0.664 


16.93 


16.80 


(1) 


C11322+3029 


0.697 


16.90 


16.72 


(1) 


C10020+0407 


0.698 


17.06 


17.01 


(1) 


C11322+3027 


0.751 


16.97 


16.79 


(1) 


C12155+0334 


0.820 


17.22 


17.12 


(1) 


C11603+4313 


0.895 


17.78 


17.64 


(1) 


C11603+4329 


0.920 


18.24 


18.02 


(1) 



a Coma cluster. 

b Also known as F1557.19TC (Couch et al. 1991). 
c Also known as J1888.16CL (Couch et al. 1991). 
d Also known as F1767.10TC (Couch et al. 1991). 

References. — (1) Aragon-Salamanca et al. 1993; (2) Barger et al. 
1996; (3) Barger et al. 1997 



to the three lowest redshift points. Since we do not have K 
magnitudes for many low redshift galaxies, the normalisa- 
tion of the no-evolution line could be somewhat uncertain. 
Howev er, extensive optical photometry exists and, given the 



shift clusters from the ground-based images. Many of inter- 
mediate and high redshift clusters have been imaged with 
HST, thus morpho logical informatio n exists for 16 of the 
z > 0.3 clusters ((|Smail et al. 1996[), 1997; flCouch et aT 



Lauer 1 995[ )), we feel justified to use the i?-band photometry 
of Lauer & Postman (1992) to estimate 7f-band magnitudes. 
The shaded area at low-redshift in Figure |l| represents the 
if-band magnitudes for z < 0.06 BCGs estimated from the 
i?-band data of Postman & Lauer, corrected to our met- 
ric aperture (using d log L/d log r = 0.5; from the sa me au- 
thors), and assuming a typical R — K — 2.6 (cf. (Bower 



homogeneity in the colours of these galaxies ((Postman & 1997)). In general, the BCGs are either cD galaxies or giant 

ellipticals, although for the highest- z clusters there could be 
some ambiguity since the extended cD halos might not be 
clearly visible given the relatively high surface brightness 
limits achievable with HST. Thus some of the BCGs clas- 
sified as ellipticals could be cD galaxies. The exception to 
this rule is the BCG in the 3C295 cluster, which is a radio 
galaxy and shows an un-resolved AGN-like morphology in a 
disturbed disk or envelope. 

The optical and optical-infrared colours (from ground- 
based and HST data) of the BCGs are compatible with those 
of the colour-magnitude sequence of cluster ellipticals and 
SOs. The exception is, again, 3C295, which is substantially 
bluer (by ~ 0.8 mag in R — K). The peculiar colours and 
morphology of this galaxy are probably related to its being 
a powerful radio source. However, its Jf-band luminosity 
agrees very well with that of the rest of the BCGs in our 
sample, and we have kept it in our analysis. Taking it out 
would not alter our conclusions. 



Lucey , l Ellis 1992a ),b). The width of the shaded area re- 
flects the observed scatter. The agreement with our local 
normalisation is remarkably good. 

In agreement with the results of Aragon-Salamanca 
et al., the scatter of the observed magnitudes around the 
no-evolution prediction is 0.30 mag. Assuming that the K- 
band light provides an estimate of the total stellar mass 
of the galaxies, this very low scatter implies that BCGs at 
a given redshift have very similar masses in stars (within 
30% r.m.s.), which has yet to be explained by any model of 
galaxy formation (see section 3). 

Morphological information is available for the low red- 
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Figure 1. (a) Magnitude-redshift relation (Hubble diagram) for 
the Brightest Cluster Galaxies in the rest-frame i^-band. The 
K magnitudes have been measured inside a projected 50kpc- 
diameter aperture. The solid line shows the no-evolution predic- 
tion, normalised to the low-redshift data. The dashed line corre- 
sponds to the luminosity evolution expected for a stellar popula- 
tion formed at 2f or = 2 evolving passively. The dotted line shows 
the evolution expected for Zf OI = 5. The shaded area at low- 
redshift represents the i\-band magnitudes of z < 0.06 BCGs 
estimated from the i?,-band data of Postman & Lauer. The width 
of the shaded area reflects the observed scatter (see text for de- 
tails). Ho = 50kms — 1 Mpc _1 and go = 0.0 were assumed, (b) 
As (a) but for go = 0.5. 

2.2 Interpretation using evolutionary synthesis 
models 

Figure [j] shows that the luminosity of the BCGs (inside a 
fixed metric aperture of 50 kpc diameter) does not evolve 
strongly with redshift. This is shown more clearly on the 
top panel of Figure ^| where the no-evolution prediction has 
been subtracted from the data. It is clear that for qo = 0.0 
the observed magnitudes do not seem to evolve with redshift. 
For qo = 0.5 there is a hint of negative evolution: BCGs at 
high redshift tend to be marginally fainter than low redshift 
ones. 



Figure 2. (a) Top panel: The same data presented in figure |l| 
after subtracting the no-evolution prediction. Middle and bottom 
panels: the same data after subtracting models for luminosity 
evolution in which the BCG stars form at Z{ OI = 2 and Zf or = 5 
respectively and evolve passively thereafter (for go = 0.0). The 
solid lines represent least-squares linear fits, (b) As (a) but for 
go = 0.5. 

A convenient parametrisation of the luminosity evolu- 
tion is Lk(z) = Lk{0) x (1 + z) 7 . Least-squares fits to the 
data yield 7 = -0.06 ± 0.3 and -0.6 ± 0.3, for qo = 0.0 and 
0.5 respectively. Thus the BCGs show no or marginally neg- 
ative luminosity evolution. However, as mentioned above, 
the colours of these galaxies show the same evolution as the 
early-type cluster galaxies: they get progressively bluer with 
redshift at a rate which indicates that their stellar popula- 
tions were formed at z > 2 and evolved passively thereafter. 
If the total stellar mass of the galaxies has remained con- 
stant, we would expect them to get progressively brighter 
with redshift, as the average ages of their stellar popula- 
tions get younger. Since that brightening is not observed, 
the most likely explanation is that the total mass in stars 
in the BCGs has not remained constant but has grown with 
time. We will now estimate the rate of this growth. 

We used evolutionary population synthesis models 
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(( Bruzual fc Chariot 1997 ); see also (Bruzual & Chariot 
1993 |) ) | to predict the expected luminosity evolution of a 
passively- evolving stellar population formed at a given red- 
shift. That should represent the brightening of the stellar 
populations due to the decrease in average stellar age with 
Spain CI OS 



redshift) . A Scalo (1086) Initial Maoo Function with atcllar 



masses l in the range O.IMq < m < 125 A/® was assumed, 

The exact value of the upper mass limit is not critical given 
the expected ages of the stellar populations in the observed 
redshiftr 



range, provided that it ia above the turn off maoa 
at the AgoB corresponding to z < 1, The model predictions 



are plotted on Figure [j] (LE lines) for formation redshifts 
Zfor = 2 and Z[ OT = 5. These models produce a colour evolu- 
tion compatible with the observations of early-type cluster 
galaxies. Note that Zf or < 2 is clearly ruled out by the ob- 
served slow colour evolution (cf. section 1). Making zf OI > 5 
makes very little difference to the age of the stars. We have 
assumed solar metallicity, but changing the metal content 
does not significantly alter the predictions. Other mode ls 



(e.g. (Arimoto fc Yoshii 1986; Kodama fc Arimoto 1997)) 



predict very similar luminosity evolution, thus we believe 
that the models considered here bracket reasonably well the 
luminosity evolution expected from the colour constraints 
and that the predictions are not very model-dependent. 

In Figure ^ (lower two panels) we have plotted the ob- 
served magnitudes after subtracting the luminosity evolu- 
tion predictions. Since the models assume a constant stellar 
mass, the rate of change with redshift in the K magnitudes, 
after correcting for the expected brightening due to passive 
evolution, should be a direct measurement of the rate of 
change in stellar mass. Note that strictly speaking, what we 
call "mass in stars" really means "mass in luminous stars" , 
i.e., those which contribute to the galaxy A"-band light. Pa- 
ramete tiEing tho ovolution ae M*(z) — M*(0) x (1 — a) 7 , 



least- gqliiarpg fits ghrp trip fnllnwing rpgnltg- 



go =0.0 zfor = 2 7 = -1.2 ±0.3 

Zfor = 5 7 = -0.7 ±0.3 

g = 0.5 ^or = 2 7 = -2.2 ± 0.3 

Zfor = 5 7 = -1.7 ±0.3 

These rates imply that the mass in stars of a typical 
BCG grew by a factor ~ 2 if go = 0.0 or ~ 4 if go = 0.5 from 
z ~ 1 to z ~ 0, i.e., in the last ~ 8-10 Gyrs. 

At this point we would like to make a historical com- 
ment. One of us (AAS) carried out the analysis described in 
this section independently of the rest of the authors and in- 
vited the Durham and Munich groups to make predictions, 
using the semi-analytic models described in section 3, for 
the growth of the stellar mass in BCGs as a function of 
redshift. In the first instance, this was done without any 
prior knowledge of the observationally determined rates and 
without any interaction between the two groups to ensure 
that the predictions where truly independent. The surpris- 
ing level of agreement of the model predictions with each 
other, and with the observations, encouraged us to compare 
the model output with the data in more detail. This is pre- 
sented in the next section. 



3 GALAXY FORMATION AND EVOLUTION 
MODELS 

In this section, we discuss the predictions of the semi- 
analytic models of the Durham and Munich g roups. Fu ll 



detai ls of t he Durham models can be f ound in (Cole et al 



1994) and (Baugh, Cole fc Frenk 19961). The version of the 



Munich models used here is the same as that described in 



(Kauffmann & Chariot 1997). Further details about the Mu- 



nich mode ls can be found in (Kauffmann, White & Guider- 
doni 1993j ). Although the two models are very similar in 



outline and in their basic framework, many features, for ex- 
ample the simple parameterisations used to describe star 
formation and feedback, are different in detail. 

Briefly, the hierarchical collapse and merging of an en- 
semble of dark matter halos is followed using Monte Carlo 
techniques. Gas virialises in the dark matter halos, cools ra- 
diatively and condenses onto a central galaxy at the core of 
each halo. Star formation occurs at a rate proportional to 
the mass of cold gas present. The supply of cold gas is regu- 
lated by feedback processes associated with star formation, 
such as stellar winds and supernovae. As time proceeds, a 
halo will merge with a number of others, forming a new halo 
of larger mass. All gas which has not already cooled is as- 
sumed to be shock heated to the virial temperature of this 
new halo. This hot gas then cools onto the central galaxy of 
the new halo, which is identified with the central galaxy of 
its largest progenitor. The central galaxies of the other pro- 
genitors become satellite galaxies, which are able to merge 
with the central galaxy on a dynamical friction timescale. If 
a merger takes place between two galaxies of roughly compa- 
rable mass, the merger remnant is labelled as an "elliptical" 
and all cold gas is transformed instantaneously into stars in 
a "starburst". T he same stellar population models ((Bruzual 



Chariot 1997)) used in Section 2.2 are employed to turn 



the star formation histories of the model galaxies into broad 
band luminosities. 

In these models, there are three different ways in which 
the stellar masses of the brightest galaxies in clusters can 
grow: 

(i) merging of satellite galaxies that sink to the center of 
the halo through dynamical friction 

(ii) quiescent star formation as a result of gas cooling from 
the surrounding hot halo medium 

(iii) "bursts" of star formation associated with the accre- 
tion of a massive satellite. 



As noted by (Kauffmann, White & Guiderdoni 1993) 



the colours and absolute magnitudes of central cluster galax- 
ies are predicted to be bluer and brighter than observed if 
all the gas present in the cooling flows of massive clusters 
turns into visible stars. To fix this problem, these authors 
simply switched off star formation in cooling flows when the 
circular velocity of the halo exceeded 500 kms -1 . This value 
was chosen so as to obtain a good fit to the bright end of 
the Virgo cluster luminosity function. Another solution is 
to assume that gas in the central regions of halos follows a 
different density profile to that of the dark matter. In partic- 
ular, if the gas has a constant density core, cooling is very 
much less efficient in the centres of massive clusters. The 
precise form of the gas density profile is much less impor- 
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Figure 3. The build up in stellar mass of the largest progenitor 
in four examples of BCGs taken from the Munich models, which 
reside in dark matter halos with circular velocity lOOOkms - 1 at 
the present day. Star formation is switched off in any progenitor 
halo whose circular velocity is in excess of SOOkms -1 . The solid 
lines show the accumulation of stellar mass from merging events, 
the dashed lines show the mass contributed by stars forming from 
gas cooling from the surrounding hot halo medium, and the dot- 
ted lines show the mass contributed by star formation bursts as- 
sociated with major mergers. 



tant for cooling in low mass halos, since these objects are 
denser and have much shorter cooling times. 

In figure ^, we consider the evolution of a subset of the 
z — BCGs in a model with star formation switched off 
in halos with circular velocity greater than 500 km s - . The 
solid line shows the cumulative growth in stellar mass of 
the BCG due to merging events. The dashed line shows the 
cumulative contribution from stars formed from quiescently 
cooling gas, and the dotted line is the contribution from stars 
formed in bursts during major mergers. As can be seen, star 



formati an from bursts and cooling gas account for only a few 



percent of the final mass of the BCG. The rest comes from 
accreted galaxies. The four panels in the figure represent 
independent Monte Carlo realizations of the formation of a 
BCG in a halo of circular velocity lOOOkms -1 . It is striking 
that although each BCG has a significantly different merging 
history, their final stellar masses differ very little from one 
realization to another. 

We will now make a comparison with the data presented 
in section 2. The observational selection picks out the rich- 
est clusters at each redshift, but the sample is not complete 
in any well-defined sense. We mimic the observational selec- 



Figure 4. The average stellar mass of BCGs as a function of 
redshift. We have selected galaxies at each epoch that reside in 
halos with circular velocities greater than 700 km s" 1 . The points 
have been normalised by dividing by the average stellar mass of 
the BCGs in each model at z = 0. (a) shows the predictions for 
a Cold Dark Matter universe with the critical density, whilst (b) 
shows an open universe with a present day density parameter 
of f2o = 0.2. In each panel, the filled circles show the results 
from the Munich models and the open traingles show the results 
of the Durham models. The curves show the trends in stellar 
mass deduced in Section 2.2; the dotted curve corresponds to 
z for = 5 an d the dashed curve to Zf or = 2. The open circles in 
(a) and (b) show the evolution in the average stellar mass of the 
BCGs predicted by the Munich model when a cut in halo mass 
of 2 X 10 14 h~ 1 Mq is applied to select clusters at each redshift. 



objects with increasing redshift. We also explore the effect 
of selecting clusters by mass instead of circular velocity. 

We present the model predictions for two Cold Dark 
Matter cosmologies; one with the critical density and an 
open model with density parameter Ho = 0.2. The density 
fluctuations in each case are n ormalised to reproduce the 



abundance of rich clusters (e.g. (White, Efstathiou fc Frcnk 



1993); (Eke, Cole fc Frenk, 1996); (Viana fc Liddle 1996)) 



At several redshifts in the range < z < 1, we have selected 
~ 200 halos with circular velocities greater than 700 kms -1 . 
Note that the clusters are weighted so that they are rep- 
resentative of the mass distribution of halos at the given 
redshift. 

The evolution of the average stellar mass of the bright- 
est galaxies in the halos with redshift is shown in Figure ^j. 
We have divided the average mass of the BCGs at each red- 
shift by the average mass of BCGs found in the redshift zero 
clusters. The open triangles show the results of the Durham 



models, which adopt the halo density profile of (Navarro 



tion as | best we can by calculating the masses of BCGs in Frenk fc White 199q ) for the dark matter. The gas is as- 



halos with circular velocities greater than some fixed value 
as a function of redshift. This means we are selecting rarer 



sumed to follow a different density profile, with a constant 
density at the core of the halo (full details of these extensions 
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Figure 5. The rest frame V — I colour of the model BCGs as 
a function of redshift. Again, (a) shows the results for a critical 
density CDM universe, whilst (b) shows the model output for an 
open universe with Qq = 0.2. The filled symbols show the Munich 
models and the open symbols show the Durham models. The 
open stars in (a) show a Durham model in which the gas density 
profile is the same as that of the dark matter in the halo. Hence, 
gas cools more efficiently and more recent star formation takes 
place, making the BCG's bluer. The line is a fit to the observed 
colour evolution of the BCG's taken from Aragon-Salamanca et 
al. (1993). 
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Figure 6. (a) The average stellar mass of BCGs at z = as a 
function of the circular velocity cut used to select clusters. The 
mean BCG mass increases by a factor of two over a range in 
circular velocity that corresponds to a increase by a factor of 5 
in dark matter halo mass. The filled symbols show the Munich 
results and the open symbols show the Durham results, (b) The 
scatter in the absolute K-magnitudes of BCGs as a function of the 
circular velocity used to define the cluster sample. The Durham 
models predict a larger scatter than the Munich models. 



to the original flCole et al. 1994[ ) model will be given in Cole 
et al, in preparation). The filled circles show the predictions 
of the Munich models, in which isothermal density profiles 
are assumed for both the gas and the dark matter. Star for- 
mation is switched off in halos with a circular velocity in 
excess of 500kms _1 . The curves show the growth in stellar 
mass with redshift of the BCGs deduced in Section 2.2; the 
dotted lines show the parameterisation with Zf or — 5 whilst 
the dashed line shows Zf or — 2. Both models fit the observed 
evolution very well. 

In Figure [| we show the evolution in the average rest- 
frame V — I colours of the BCGs. Once again, the open 
triangles show the Durham models and the filled circles the 
Munich models. The curves are a fit to the observed colour 



evolution given in (Aragon-Salamanca et al. 1993). For the 



S7 = 1 cosmology, both models fit the observed colour evo- 
lution reasonably well. The Durham models predict larger 
scatter in the colours because the BCGs are still undergo- 
ing some low-level star formation, whereas the BCGs in the 
Munich models have ceased forming stars altogether. For 
reference, we also show a model (shown by the open stars) 
where cooling and star formation are efficient at the cen- 
tres of massive clusters. This is a model where both the gas 
and th e dar k matter follow the standard (Navarro, Frenk & 



996) density profile. In this case, the colours are too 
0.2 magnitudes. In the low-density cosmol- 



White 

blue by ~ 0.1 

ogy, it is interesting that the Durham and Munich models 



exhibit opposite trends in colour. The Munich models are 
redder than in a flat cosmology, because the BCGs form ear- 
lier and contain older stars. On the other hand, the Durham 
models are bluer. This is because halos in a low-density Uni- 
verse are more concentrated and gas still cools at the cluster 
centres, even if a different density profile is assumed. 

An important question to consider is how our results 
would be affected by adopting a different selection criterion 
for the clusters in the model. To explore this, Figure 
shows how the average mass of a z = BCG changes as 
a function of circular velocity cutoff for both the Munich 
and Durham models. This plot shows that more massive 
BCGs are found in more massive halos. The stellar mass of 
a BCG changes by a factor of 2 for a factor 5 increase in halo 
mass. Figure shows how the predicted scatter in K-band 
absolute magnitude changes as a function of circular velocity 
cutoff. As can be seen, the scatter decreases as richer clusters 
are selected. Both the Munich and Durham models predict 
larger scatter than that observed, though the discrepancy is 
worse for the Durham models. 

It should be noted that clusters of fixed circular velocity 



are less massive at high redshift by a factor : 



N3/2 



(this 



scaling arises from the fact that halos collapse and virialize 
at a density ~ 200p(z), where p(z) is the mean density of 
the Universe at redshift z). Selecting clusters according to a 
halo mass cutoff, rather than a circular velocity cutoff, will 
thus select more massive BCGs at high redshift. The effect 
of mass selection on the predicted evolution is shown by 
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the open circles in Figure ^, where we have selected clusters 
greater than 2x 1O 14 M0 at each redshift. The mass evolution 
is now slightly weaker, but the effect is small and the model 
results are still in good agreement with the observations. 



4 DISCUSSION 

We have shown that the 7f-band Hubble diagram for a 
sample of brightest cluster galaxies in the redshift range 
< z < 1 has a very small scatter (0.3 magnitudes), in good 
agreement with earlier studies. The BCGs exhibit very little 
luminosity evolution in this redshift range: if go = 0.0 we de- 
tect no luminosity evolution; for go = 0.5 we measure a small 
negative evolution (i.e., BCGs were about 0.5 magnitudes 
fainter at z = 1 than today). But substantial positive lu- 
minosity evolution would be expected if the mass in stars 
of these galaxies had remained constant over this period of 
time: BCGs should have been brighter in the past since their 
stars were younger. The most likely explanation for the ob- 
served zero or negative evolution is that the stellar mass of 
the BCGs has been assembled over time through merging 
and accretion, as expected in hierarchical models of galaxy 
formation. Since the colour evolution of the BCGs is con- 
sistent with that of an old stellar population (zf orrn > 2) 
that is evolving passively, we have used evolutionary pop- 
ulation synthesis models to estimate the rate of growth in 
stellar mass for these systems. We find that the stellar mass 
in a typical BCG has grown by a factor ~ 2 since z ~ 1 if 
go = 0.0 or by factor ~ 4 if go = 0.5. 

Two independent groups of galaxy formation modellers 
predicted the decrease in stellar mass of the BCGs with in- 
creasing redshift, in advance of seeing the data and without 
any fine tuning of their semi-analytic model parameters. We 
have extracted the properties of BCGs in massive halos with 
circular velocities greater than 700 km s _1 or masses greater 
than 2 x 10 14 M Q , at various redshifts. The average stellar 
mass of the BCGs increases by a factor of ~ 4 — 5 in the crit- 
ical density models. In the low density models, the growth 
in stellar mass is more modest, between a factor 2 to 3. This 
is because clusters assemble at higher redshift in a low den- 
sity Universe. Note, however, that both models agree well 
with the data when the appropriate value of go is used in 
the analysis. It is not possible to make statements about a 
preferred cosmology from this exercise. 

We also find that in order to reproduce the observed 
colours of the galaxies, some mechanism must act to sup- 
press visible star formation in the cooling flows of the most 
massive halos. It has been postulated that much of this 
gas may be for ming low-mass stars U FjdDkm i ^Nulsen& 
Canizares, 1982] )) or cold gas clouds ((White et al. 1991)). 
Another possibility we have explored is that the gas may 
be less concentrated in the centres of clusters than the dark 
matter, and hence not able to cool efficiently. In both cases, 
the evolution in the _K"-band luminosities of the BCGs comes 
about almost exclusively as a result of the merging of satel- 
lite galaxies in the cluster. An important and unanticipated 
result of our analysis is that the predicted scatter in the 
stellar masses of the BCGs built through merging and ac- 
cretion is still quite small, although somewhat larger than 
the observed one. The scatter depends on the criteria used 
to select the model clusters. 
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